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Abstract  -  The  Radar  Cross  Section  of  missiles  at  HF  band  is 
affected  by  the  presence  of  the  missile  plume.  In  fact,  although 
the  plume  is  transparent  at  the  microwave  frequencies,  it 
reflects  almost  the  totality  of  the  energy  at  the  HF  band. 

In  this  paper  we  evaluate  the  effect  of  the  missile  plume  on  its 
Radar  Cross  Section  (RCS)  at  HF  band.  In  order  to  achieve  the 
goal  we:  1)  verify  that  the  missile  plume  can  be  physically 
represented  by  a  plasma  through  a  chemical  analysis  of  the 
combustion;  2)  define  the  time -varying  non-linear  differential 
equation  that  rules  the  electron  density  equilibrium  in  the 
plume;  3)  solve  the  time -varying  non-linear  differential  equation 
in  order  to  evaluate  the  length  of  the  plume  that  contributes  to 
the  RCS  calculation;  4)  evaluate  the  RCS  of  the  missile+plume 
by  means  of  the  Method  of  Moments  (MoM).  As  a  numerical 
case,  the  application  of  both  a  surface  and  sky  wave  Over  The 
Horizon  (OTH)  radar  has  been  considered.  Specifically,  the 
RCS  (as  a  function  of  the  incidence  angle)  of  the  missile  during 
the  boost  phase  trajectory  has  been  evaluated. 

I.Introduction 

Detection  of  ballistic  missiles  in  the  boost  phase  by  Over  the 
Horizon  (OTH)  HF  radar  is  a  very  important  task  in  a  war 
scenarios  to  early  activate  the  defence  system.  The  visibility 
of  a  missile  at  long  distance  is  mainly  affected  by  its  RCS, 
which  is  a  parameter  extremely  important  to  correctly  design 
the  radar  system.  At  FF  band,  the  calculation  of  the  RCS 
obtained  by  considering  only  the  body  of  the  missile  would 
lead  to  an  underestimated  value.  In  fact,  for  frequencies 
ranging  from  3  to  30  MHz  and  considering  the  boost  phase, 
the  plume  behaves  as  plasma.  To  demonstrate  that  the  plasma 
model  is  appropriate,  the  chemical  reactions  of  the 
propulsion  system  combustion  must  be  considered1. 
Specifically,  the  exothermic  combustion  process  produces  a 
rapid  temperature  augment  that  causes  a  split  of  the  residual 
alkaline  metals  in  positive  and  negative  ions.  Negative  ions 
combine  with  the  hydrogen  and  produce  neutral  molecules 
and  electrons.  At  the  end  of  the  process,  the  exhaust  gas  is 
composed  of  neutral  molecules  and  of  a  cloud  of 
macroscopically  neutral  (“ Almost  neutrality  condition ”) 
positive  and  negative  charges.  This  gas  status  can  be  defined 
as  plasma,  confirming  the  hypothesis. 


1  In  this  paper  solid  fuel  rockets  are  considered .  For  liquid  fuel  rockets 
different  reactions  should  be  considered. 


As  already  known  in  the  case  of  the  ionosphere  effect  on  the 
radio  waves,  a  reflection  of  the  incident  e.m.  field  occurs  at 
suitable  frequencies  and  incidence  angles.  Therefore,  the 
radar  echo  energy  is  in  general  affected  by  the  presence  of 
the  missile  plume. 

The  aim  of  the  paper  is  to  evaluate  the  contribution  of  the 
missile  plume  to  its  RCS.  To  obtain  this  result  we  1)  define 
the  time-varying  non-linear  differential  equation  that  rules 
the  electron  density  equilibrium  in  the  plume;  2)  numerically 
solve  the  differential  equation;  3)  apply  the  Method  of 
Moments  (MoM)  to  the  composition  of  missile  body  and 
plume,  both  modeled  as  cylinders.  The  numerical  calculation 
of  the  RCS  of  a  short  range  ballistic  missile  (length  of  6  m) 
at  4  different  frequencies  and  with  incident  angles  (see  Oj  in 
fig.  1)  ranging  from  90°  to  0°  is  presented  and  discussed  in 
the  last  section. 

n.  Plume  electromagnetic  model 

In  first  part  of  this  section  the  chemical  reactions  due  to  the 
gas  post-combustion  are  studied.  By  following  chemical 
principles,  the  calculation  of  the  electron  concentration  at 
the  missile  nozzle  is  then  evaluated.  In  the  second  part  of  the 
section,  the  time-varying  non-linear  differential  equation  that 
rules  the  electron  concentration  in  the  missile  plume  is  first 
defined  and  then  solved. 

II.  1  Combustion  chemical  reactions 

Various  mechanisms  occur  in  the  post-combustion  chemical 
reactions  [1]  (see  Table  I).  The  most  important  phenomena 
affecting  the  electron  density  are  the  charge  production  and 
charge  recombination.  In  this  section  we  focus  on  the 
chemical  reactions  of  a  solid  fuel  missile.  For  missile 
altitudes  lower  than  20  Km  with  respect  to  the  sea  level,  the 
main  electronic  production  mechanism  is  showed  in  eq.  1 : 

K  +  C1^K+  +  CT  (1) 

followed  by: 

H  +  Cl  ^  HC1  +  e"  (2) 

The  dominant  recombination  reaction  is  showed  in  eq.  3 : 

K+  +  e_  +  M  ^  K  +  M  (3) 

It  is  worth  noting  that,  although  the  magnesium  particles 
involved  in  the  reaction  of  eq.  3  (originally  present  in  the 
fuel)  seem  to  not  affect  the  production  of  electrons,  they 
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must  be  accounted  for  in  order  to  satisfy  impulse  and  energy 
conservation  law. 

The  reactions  showed  in  eq.  1-3  allow  us  to  calculate  the 
electron  density  at  the  missile  nozzle,  namely  n,. 
Specifically,  we  consider  the  molar  fraction  of  chlorine  in 
the  combustion  products,  which  is  typically  equal  to 
A/c/=2.7*10"3.  Therefore,  a  mole  of  fuel  contains 
Ne=Nci*Nav=  16.26  *1020  electrons/mole,  where  Nav  is  the 
Avogadro  number. 


Reactions 

Chemical 

components 

Products 

Hydrogen  Carbon 

O+O+M 

o2+m 

Oxide  combustion 

O+H+M 

OH+M 

H+H+M 

h2+m 

H+OH+M 

h2o+m 

CO+O+M 

co2+m 

OH+OH 

h2o+o 

OH+  H2 

h2o+h 

0+  h2 

OH+H 

H+  02 

OH+H 

CO+OH 

co2+h 

Reactions  containing 

H+Cl+M 

HCl+M 

chlorine 

C1+  h2 

HC1+H 

h2o+ci 

HC1+OH 

OH+C1 

HCl+O 

Neutral  Potassium 

K+HC1 

KC1+H 

Electrons  production 

K++  e‘+M 

K+M 

and  recombination 

K++cr 

K+Cl 

Charge  exchange 

C1+M+  e- 

M+Cl 

HC1+  e" 

H+  Cl 

Table  I  -  Main  post-combustion  chemical  reactions 


To  evaluate  ni9  the  volume  of  a  mole  of  exhaust  gas  at  the 
pressure  and  temperature  present  at  the  burner  exit  must  be 
evaluated.  Hence,  in  order  to  calculate  the  pressure  at  the 
burner  exit,  we  should  evaluate  the  pressure  distribution 
along  the  plume  by  assuming  an  adiabatic  expansion  (this 
hypothesis  is  realistic  because  the  pressure  gradient  should 
be  intuitively  higher  than  the  temperature  gradient).  Although 
this  approach  is  physically  consistent,  we  prefer  to 
drastically  simplify  the  problem  by  assuming  that  the 
electron  density  is  not  affected  by  variations  in  fluid 
pressure.  Therefore,  we  make  the  assumption  that  it  only 
depends  on  electron-ion  recombination.  In  this  condition,  the 
electron  density  at  the  burner  exit  is  assumed  to  be  equal  to 
the  one  at  great  distance,  where  the  recombination  effect  is 
negligible.  In  other  words,  we  evaluate  rii  assuming  the 
pressure  at  the  burner  exit  equal  to  the  atmospheric  pressure. 
Although  this  assumption  produces  and  underestimation  of 
rii,  it  does  not  significantly  affect  the  evaluation  of  the 
missile  RCS.  This  statement  will  be  demonstrated  in  section 

m. 

Assuming  that  the  burner  exit  temperature  is  equal  to  2000  K 
[1]  and  applying  the  perfect  gas  equation,  we  obtain  a  value  of 


the  volume  V  =  nRT/P  =  16.629*10  m  ,  where  R  is  the 
Boltzmann  constant  and  n  is  the  number  of  moles  (set  to  1). 
In  formula,  the  electron  density  nt  is  showed  in  eq.  4: 

N  u  i 

n.  =  — —  =  10  cm  3  (4) 

*  y 

It  is  worth  noting  that  when  a  solid  fuel  is  used  the  chemical 
reactions  change.  Nevertheless,  the  electron  density  at  the 
burner  exit  can  be  evaluated  by  means  of  an  equivalent 
chemical  analysis. 

II.2  Plume  electron  density  evaluation 

Referring  to  Fig.l,  let  the  origin  of  a  reference  system  be 
located  at  the  center  of  the  missile  burner  nozzle,  and  the  x 
coordinate  be  directed  along  the  missile  axis. 


Fig.l  -  System  geometry 


If  we  assume  that  for  each  value  of  x,  the  electron  density 
is  uniformly  distributed  on  a  circular  plane  with  diameter  D 
equal  to  the  one  of  the  missile2,  the  differential  equation  that 
rules  the  spatial  and  time  evolution  of  the  electron  density 
ne  ( t ,  x)  is  given  by  eq.  5 : 

d  n  (t,x)  9 

— £ - h  ocn2e(t,x)  =  q(t,  x)  (5) 

dt 

where  q(t,  x)  is  the  forcing  term  for  the  generic  section  x,  v 

is  the  electron  conveying  speed  and  a  is  the  recombination 
factor.  Note  that  the  eq.  5,  when  the  forcing  term  is  equal  to 
zero  ( q(t,  x)  =  0),  represents  the  electron  recombination  law 
for  a  generic  plasma  [2].  The  definition  of  the  forcing  term 
will  be  further  given.  It  is  worth  noting  that  the  electron 
density  ne(t,x)  is  non-zero  after  T=x/v  seconds.  In  fact,  T 

represents  the  time  interval  necessary  for  the  electrons  to 
reach  the  coordinate  x  from  the  origin.  For  the  sake  of 
simplicity,  let  the  coordinate  x  be  divided  in  M  sub-cells:  let 
be  xm  =  m  Ax  where  m= 0, 1 ,2, . . . ,  M- 1 .  The  quantity  L-MAx 
represents  the  plume  length  to  be  analyzed.  By  denoting  with 
ne{t’xm)  =  ne,m  (0  the  electron  density  of  the  ra-th  sub-cell, 


2 

As  a  first  approximation  it  is  possible  to  neglect  the  quantity  of  electrons 
that  exit  the  cylinder. 
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the  differential  equation  that  rules  the  equilibrium  in  the  m-th 
sub-cell  can  be  written  as  follows: 
d  n  (t)  0 

+  an2em(t)  =  q(t,xm )  (6) 

at 

with  the  condition:  ne  m(t)  =  0  when  t<Tm  -  xm/v .  To 

completely  describe  eq.  6  we  have  to  define  the  forcing  term 
q(t,xm)  and  the  recombination  coefficient  a. 

The  forcing  term  q{t,xm )  accounts  for  the  variation  of  the 
electron  concentration  due  to  the  electrons  propagation 
along  the  x  axis.  In  fact,  if  we  apply  the  mass  conservation 
principle  to  the  generic  m-th  sub-cell,  we  have  to  consider 
the  nunmber  of  electrons  that  enter  and  exit  the  sub- cell.  In 
order  to  calculate  the  variation  of  the  electron  concentration 
per  time  unit  we  refer  to  fig.  2. 


Fig.2  -  Sub -cell  electron  density  equilibrium 


We  first  consider  the  interface  between  the  m-th  and  (m+1)- 
th  sub-cells.  Let  SV  be  the  infinitive  volume  of  a  cylinder  of 
diameter  D  and  height  8x  that  includes  the  interface  between 
the  two  sub- cells.  The  infinitive  height  8x  can  be  written  as 
Sx=v  St ,  where  St  is  an  infinitive  time  interval  that  allows  us 
to  consider  the  electron  concentration  time -unvarying.  The 
volume  8V  can  be  calculated  through  eq.  7 : 

7i  D2 

8v  =  -——vSt  (7) 

4 

The  number  of  electrons  that  exit  the  m-th  sub-cell  (enter 
the  (m+l)-th  sub-cell)  in  the  infinitive  time  interval  St  can  be 
written  as  follows: 


(8) 

In  order  to  consider  the  variation  of  the  electron 
concentration  in  the  sub-cell,  the  number  of  electrons  that 
exit  the  m-th  sub- cell  must  be: 

1)  divided  by  the  volume  of  the  sub-cell,  which  is  equal  to 


V: 


kDl 


Ax ; 


2)  divided  by  the  infinitive  time  interval  St  and  take  the 
limit  for  &-> 0,  as  follows: 


d  r 

—  n  ^  m 

dt 1  e’m^w+1 


(t)~}  =  lim  — 1 — N  ( t)  =  —n  (t) 

v  J  &-*>  St  V  e’m  w  Ax  e’m  V  ' 


(9) 


The  calculation  of  the  number  of  electrons  that  enter  the  m- 
th  sub-cell  (exit  the  (m-l)-th  sub-cell)  is  equivalent.  The  only 
exception  is  represented  by  the  concentration  of  the 
electrons,  which  is  evaluated  for  the  (m-l)-th  sub-cell  at  the 
time  t. 

By  avoiding  mathematical  passages,  we  can  write: 

Yt  (  *)]  =  ^  “e.m-l  (0  (  1  0) 

Hence,  the  forcing  term  q(t,xm)  can  be  written  as  follows: 

q{^xm)= (*)-"«,<»  M]  C1 !) 

By  combining  eq.  6  and  eq.  11  we  obtain  the  differential 
equation  that  rules  the  equilibrium  of  the  electron 
concentration  for  the  generic  m-th  sub-cell,  as  follows  in  eq. 
12: 


d  ne  m  it) 
dt 


2  V 

(Xtl  m  (0  - ne 

Ax  ’ 


,(o=L^-1(o 

Ax 


(12) 


with  the  condition:  ne  m(t)  =  0  when  t  <  Tm  =  xm/v  . 


The  solution  of  the  non-linear  differential  equation  with 
time- varying  coefficients  of  eq.  12  can  be  obtained  through 
an  iterative  numerical  algorithm.  Specifically,  considering 
the  first  sub-cell  (m=0),  we  have  to  solve  the  following 
differential  equation: 


dne  0(0 

dt 


+  CC  Yle  q  {t)  + 


Ax 


,(0  = 


Ax 


(13) 


with  the  condition:  ne  0(t)  =  0  when  t<  0.  The  initial 

concentration  rii  is  given  in  eq.  4.  In  the  second  sub-cell 
(m=l),  the  differential  equation  to  be  solved  is  the 
following: 

dn  At)  9  v  /  x  v  ,  x 
— j - +  anel(t)+— nsl(t)=—ne0(t)  (14) 


with  the  condition:  ne  l(t)  =  0  when  t  <  Tx  =  xl  /  v .  It  is  worth 

noting  that  the  term  at  the  second  member  has  been 
calculated  at  the  previous  step.  This  procedure  is  then 
iterated  for  the  next  sub-cells  (m=2,3,. .  .,M-1). 

To  solve  eq.  12,  three  parameters  must  be  defined: 

1)  The  electron  conveying  speed  v; 

2)  The  cell  length  Ax  ; 

3)  The  recombination  coefficient  a. 

The  electron  conveying  speed  v  is  a  parameter  that  depends 
on  the  missile  propulsion  system  and  it  must  be  evaluated 
case  by  case.  Concerning  the  cell  length  problem,  to  reduce 
the  errors  introduced  by  the  x  coordinate  sampling,  we  would 
be  tempted  to  take  a  value  of  Ax  approaching  to  zero. 
Nevertheless,  the  validity  of  the  differential  equations  of  eq. 
12  is  constrained  by  the  almost  neutrality  condition. 
Specifically,  the  minimum  size  of  the  cell  must  be  larger 
than  the  Debye  radius,  in  formula: 

Ax>  rD  =  6.8678.  —  (15) 

V  N 
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where  Te  is  the  temperature  and  N  is  electron  density  in  the 
specific  cell. 

The  recombination  coefficient  a  depends  on  the  electron 
temperature  and  concentration  according  to  the  following 
relationship  [2]: 


a  - 


10“ 


-n(t) 


(16) 


Eq.  16  states  that  the  recombination  effect: 

1)  increases  when  the  electronic  density  increases  (the 
larger  the  electron  density  the  higher  the  recombination) 

2)  decreases  when  the  temperature  increases. 

In  fact,  the  higher  temperature  increases  the  kinetic  energy 
of  particles,  which  overcome  the  ions  attraction.  Therefore, 
it  produces  a  reduction  of  the  ion-electron  recombination. 

By  substituting  eq.  15  in  eq.  12  we  obtain: 


d  (0 


10“ 


(t)+—n  (t)=—n  .  (t)  (17) 

’  Ax  e’",y  ’  /Ax  e'm~lK  ’ 


dt  Te  ( m ) 

where  the  function  Te  ( m )  represents  the  temperature  profile 
along  the  discrete  v  coordinate 


III.  Missile  rcs  calculation 

In  this  section  we  propose  a  numerical  example  relative  to 
the  calculation  of  the  RCS  of  a  missile  in  the  following 
scenario: 

1)  Use  of  an  OTH  radar  (at  HF  band); 

2)  Missile  in  the  boost  phase. 

In  order  to  calculate  the  RCS  we  first  calculate  the  plume 
length  that  affect  the  RCS  evaluation  and  than  we  apply  the 
MoM  to  evaluate  the  missile+plume  RCS. 

III.  1  Plume  length  evaluation 

Before  starting  with  the  plume  length  evaluation  some 
parameters  must  be  defined.  Specifically,  the  electron 
conveying  speed  within  the  plume  cylinder  v  and  the 
temperature  profile  Te(m).  The  value  of  the  electron 
conveying  speed  v  is  set  to  be  equal  to  v=100  m/s  and  the 
temperature  profile  Te{m)  is  plotted  in  fig.  3.  It  is  worth 
noting  that  the  specific  values  of  the  two  parameters  have 
been  chosen  on  a  heuristic  basis.  The  electron  density 
distribution  can  be  obtained  y  solving  eq.  17.  In  fig.  4  we 
show  the  electron  density  distribution  over  the  v  axis  for  6 
different  time  instants.  It  is  worth  noting  that  the  regime  is 
reached  after  roughly  0.30  s.  Therefore,  it  is  correct  to 
consider  the  electron  density  profile  at  £=ls  to  evaluate  the 
plume  length.  To  this  purpose  we  introduce  the  parameter  nc, 
which  represents  the  critical  electron  density.  To  simplify 
the  calculus  of  the  plume  length  we  assume  that:  the  plume 
reflects  all  the  incident  e.m.  energy  when  the  electron 


density  is  higher  than  nc  and  is  transparent  when  the  electron 
density  is  lower  than  nc? 


2  3  4 

Distance  from  the  nozzle  x  (m) 


t=0.01  sec 
t=0.015  sec 
t=0.02  sec 
t=0.025  sec 
t=0.030  sec 
t=1  sec 


Fig.  4  -  Electron  density  ne  [t,  x)  at  various  time  instants 


The  expression  of  nc  is  showed  in  eq.  18  [2]: 

■  2/T~ 


n 


-m* 


(18) 


where  m  and  q  are  the  mass  and  charge  of  the  electron  and  £o 
is  the  free  space  permittivity.  Values  of  nc  for  frequencies  in 
the  HF  Band  are  reported  in  table  n.  The  equivalent  electric 
length  /  of  the  plume  is  calculated  by  using  the  electron 
density  profile  in  a  regime  condition,  the  analytical 
expression  is  given  by  eq.  19: 

l  =  max{x:ne(tregime,x)>nc]  (19) 


3 

It  is  worth  noting  that  the  assumption  made  is  quite  strong.  It  is  the  authors’ 
intention  to  give  an  approximated  result,  which  can  be  refined  by  taking  into 
account  the  e.m.  returns  due  to  contribution  of  the  plume  tail,  where  the 
electron  density  is  below  the  threshold. 
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Frequency 

Critical  electron 
density  nc 

Plume  length 

1 

10  MHz 

1.2439*105  cm"3 

6.377  m 

15  MHz 

1 .8658*  105  cm"3 

6.342  m 

20  MHz 

2.4878*105  cm"3 

6.307  m 

25  MHz 

3.1097*105  cm"3 

6.291  m 

Table  II  -  Plume  length  at  various  frequencies 


By  zooming  in  Fig.4  around  the  intersections  between  the 
electron  density  curve  relative  to  the  time  f=ls  and  the 
horizontal  lines  which  represent  the  values  of  nc  relative  to 
the  frequencies  of  10MHz  and  25MHz,  we  note  that  (see 
fig-5): 

1)  the  plume  length  is  roughly  the  same  for  both  the 
frequencies  (hence  for  all  the  intermediate  frequencies); 

2)  the  plume  length  is  roughly  equal  to  /  =  6.5  m. 

x  10s 


Fig.5.  -  Zoom  in  of  the  electron  density  curve  of  Fig.4 


Before  evaluating  the  RCS  of  the  missile,  it  is  worth 
discussing  the  assumption  made  in  section  2.  Specifically, 
the  initial  electronic  density  of  eq.  4  is  calculated  by 
ignoring  the  real  pressure  profile  outside  the  missile  nozzle. 


This  approximation  led  to  a  value  of  nt=  1 0 1 1 .  In  order  to 
justify  such  an  assumption  let  the  initial  electron  density  be 
m=  10144.  In  other  words,  we  assume  to  have  an  error  of  three 
magnitude  order  larger  than  the  proposed  value.  In  Fig.  6  we 
report  the  graph  of  the  electron  density  profiles  for  both  the 
values  of  n,  in  regime  condition  (t=  Is).  It  is  clearly  visible 
that  when  the  plume  reaches  the  length  of  roughly  3.6  m  the 
two  curves  are  undistinguishable.  In  particular,  when  the 
length  is  equal  to  1=6.5  m,  we  can  consider  the  plume  lengths 
identical. 

III.  2  RCS  evaluation 

To  evaluate  the  RCS  of  the  missile,  the  contribution  of  the 
body  and  plum  must  be  accounted  for.  Assuming  that  1)  the 
reflectivity  of  the  body  and  plume  are  independent  of  each 
other  and  2)  the  plume  is  as  a  perfect  electromagnetic 
conductor  (PEC),  the  RCS  of  the  missile  is  given  by  eq. 
20  [3]: 

® M  ~  O Body  +  Opium  (20) 

To  numerically  evaluate  the  RCS  of  the  missile  we  have  used 
a  commercial  electromagnetic  code  [4],  [5],  namely 
Numerical  Electromagnetic  Code  (NEC).  The  NEC 
estimates  the  RCS  of  objects  by  using  the  Method  of 
Moments.  The  results  relative  to  the  frequencies  10MHz, 
15MHz,  20MHz  and  25MHz  are  reported  in  Fig.8-11, 
respectively.  The  incident  angle  is  defined  in  fig.  1. 

We  can  note: 

1)  Within  a  wide  range  of  incident  angles  the  RCS  of  the 
missile  with  plume  is  larger  than  the  RCS  of  the  only  missile 
body.  In  particular,  when  the  frequency  is  equal  to  10  MHz, 
the  presence  of  the  plume  provokes  a  RCS  gain  of  24  dB. 
Nevertheless,  in  the  other  cases  the  RCS  is  at  least  5  dB 
larger. 

2)  At  the  frequencies  of  20,25  MHz,  at  particular  values  of 
the  incidence  angles,  the  plume  produces  a  loss  in  the  RCS 
because  of  resonance  phenomena.  In  fact,  if  we  consider  the 
curve  of  the  RCS  relative  to  25MHz  (in  fig.  11),  we  note  a 
dumping  around  a  value  of  the  incidence  angle  equal  to 
ft/=70°.  Considering  that  1)  the  projection  of  the  length  of 
the  missile  body  +  plume  to  the  axis  orthogonal  to  the 
incident  e.m.  wave  is  roughly  equal  to  ltoF=  1 1 .74  m  and  2)  the 
wavelength  correspondent  to  the  frequency  /=25MHz  is 
equal  to  12,  we  can  conclude  that  the  dumping  is  due  to  the 
resonance  effect. 

rv.  Conclusions 

The  main  contribution  of  this  paper  is  due  to  the 
development  of  a  method  for  estimating  the  RCS  of  missiles, 
taking  into  account  also  the  plume  effect.  In  order  to  achieve 
the  goal,  we  have  first  verified  that  the  missile  plume  can  be 
physically  represented  by  a  plasma  model.  In  order  to  verify 


^  ftp  1014  is  the  concentration  of  the  electrons  in  a  thermonuclear  reactor, 
hence,  it  represents  a  conservative  value 
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Missile  RCS  at  20  MHz 


the  correctness  of  the  assumption  the  chemical  reactions  of 
the  post-combustion  process  have  been  considered.  In  order 
to  calculate  the  effect  of  the  plume  on  the  missiles  RCS  we 
have  first  studied  the  equilibrium  of  the  electron  density  in 
the  plasma  and  then  defined  the  non-linear  time -varying 
differential  equation  of  the  equilibrium.  Consequently,  a 
numerical  iterative  solution  of  the  differential  equation  has 
been  carried  out.  The  use  of  the  MoM  for  the  RCS 
calculation  of  the  missile+plume  system  has  finally  led  to 
the  following  conclusions: 

1)  the  effect  of  the  plume  generally  provokes  an  increase  in 
the  value  of  the  RCS; 

2)  the  resonance  effect  can  be  studied  accordingly  with  the 
particular  system  specification. 

Next  research  steps  are: 

1)  calculate  the  effects  of  exhausts  stream,  taking  into 
account  the  spatial  distributions  of  conductivity,  permittivity 
and  magnetic  permeability; 

2)  Improve  the  geometrical  model  of  the  plume,  by  taking 
into  account  for  the  non- cylindrical  shape. 


0  10  20  30  40  50  60  70  80  90 

Angle  of  incidence 

Fig. 8  -  Missile  RCS  at  10  MHz 


0  10  20  30  40  50  60  70  80  90 
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Fig.9  -  Missile  RCS  at  15  MHz 
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Fig.  10  -  Missile  RCS  at  20  MHz 
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Fig.  1 1  -  Missile  RCS  at  25  MHz 
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